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Abstract. 

We present polarimetric 5 GHz to 43 GHz VLBI observations of the BL Lacertae object PKS 0735+178, spanning March 1996 
to May 2000. Comparison with previous and later observations suggests that the overall kinematic and structural properties 
of the jet are greatly influenced by its activity. Time intervals of enhanced activity, as reported before 1993 and after 2000 
by other studies, are followed by highly superluminal motion along a rectilinear jet. In contrast the less active state in which 
we performed our observations, shows subluminal or slow superluminal jet features propagating through a twisted jet with 
two sharp bends of about 90° within the innermost three-milliarcsecond jet structure. Proper motion estimates from the data 
presented here allow us to constrain the jet viewing angle to values < 9°, and the bulk Lorentz factor to be between 2 and 4. 

Key words, galaxies: active - galaxies: jets - polarization - BL Lacertae objects: individual: PKS 0735+ 178 - radio continuum: 
galaxies - techniques: interferometric 



1. Introduction 

PKS 0735+178 (0735+178 throughout) was first identified as a 



BL Lacertae object (BLLac) by Carswell et al. J1974J . These 
authors reported absorption lines redshifted by Zabs - 0.424, 
which was later confirmed by Rector & Stocke (2001). Given 
that only absorption lines were detected, this z a ts should be 
considered as a lower limit for the redshift of 0735+178. 
Throughout this paper we assume z = z a bs - 0.424, a Hubble 
constant H - 72kms Mpc and a Friedmann-Robertson- 
Walker cosmology with Q„, = 0.3 and Q A = 0.7. Under 
these assumptions the luminosity distance of 0735+178 is d^ = 
2263 Mpc, an angular size of one milliarcsecond (mas) on the 
plane of the sky corresponds to a linear size of 5.4 pc in the 
frame of the source, and an angular motion of 1 mas/yr corre- 
sponds to a speed of 25.1 c. 

At radio wavelengths, 0735+178 is predominantly point- 
like at arcsecond resolution (Ulvestad, Johnston & Weiler 
1983 Kollgaard et al. 1992 1. Radio images made from data 
obtained with early very long baseline interferometry (VLBI) 
arrays at intermediate centimetre wavelengths have typically 
shown a compact core and a rather straight jet extending to- 
ward the north-east (e.g., Gabuzda et al. 1994 1. Multi-epoch 
VLBI observations made before 1993 allowed the identifica- 



tion of a number of superluminal jet features, with apparent 
speeds in the range 7 c to 12 c (Gabuzda et al. 1 19941 and refer- 
ences therein). 



After 1994, VLBI observations at higher frequencies pro- 
vided evidence for a curved structure within the inner jet of 
0735+178 (Kellermann et al. fTWSl Gomez et al. fT9"9"9l . In par- 
ticular, the polarimetric 22 GHz and 43 GHz images of Gomez 
et al. (1999 1, obtained in late 1996, revealed a twisted jet 
with two sharp apparent bends within the innermost 3 mas. 
Throughout the rest of the paper, we will refer to these two 
bends as the first (~ 1 mas east of the core, where the jet turns 
toward the north) and second (~ 2 mas north-east of the core, 



Send offprint requests 
iagudo@mpifr-bonn.mpg . de 



to: 



I. 



Agudo, 



e-mail: 



where the jet turns toward the east) bends. Gomez et al. ( 1999l 
also found the polarized emission following the outermost bend 
in the jet, with the magnetic vector being parallel to the local 
jet axis. Ojha et al. (2004) showed a very similar 15 GHz po- 
larization structure for a nearby epoch. Based on 8.4 GHz and 
43 GHz VLBI observations spanning from 1995 to 1998, along 
with existing previous observations, Gomez et al. (2001 1 stud- 
ied the kinematics of the components since the first VLBI ob- 
servations of 0735+ 178, made in 1979. Their results raised the 
possibility that the jet proper motions for epochs after 1994 
could be much slower than measured before 1993. 
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Comparison of 15 GHz VLBI observations in February 
1999 and 5 GHz VSOP 1 observations in January 1999, both 
probing the same jet angular scales, showed markedly differ- 
ent jet structures at these two frequencies (Gabuzda, Gomez 
& Agudo 2001 1. The most striking difference was the relative 
weakness of the 5 GHz emission near the first bend, which im- 
plied a highly inverted spectrum in this area. This behaviour 
was shown to be consistent with free-free absorption of the 
emission near this bend. Evidence for a local enhancement in 
the Faraday rotation measure (RM) of (120 + 55)rad/m 2 was 
found for the same region. This behaviour was interpreted as 
being consistent with an interaction of the jet with the external 
medium surrounding the region of the first sharp bend, similar 
to the situation observed for the VLBI jet of the radio galaxy 
3C 120 (Gomez et al. lIUUOl 

Here we present polarimetric 5 GHz to 43 GHz VLBI ob- 
servations of 0735+178, spanning from March 1996 to May 
2000. This allows us to trace the magnetic field pattern where 
the jet executes the two bends. We also carry out a more com- 
plete study of proper motions and changes in intensity along 
the jet than has been possible in earlier studies. 

2. Observations and data reduction 

The new polarimetric observations presented in this paper 
were obtained within two observing programs, both using the 
Very Long Baseline Array (VLBA) of the National Radio 
Astronomy Observatory (NRAO"). The first consisted of four 
observing sessions of 12 h performed on 26 March 1996, 2 
October 1996, 6 April 1997 and 18 October 1997 at 8.4 GHz 
and 22 GHz. During the 6 April 1997 observations, 5 GHz and 
15 GHz were also measured, see Tabled The second observ- 
ing program consisted of three observing sessions of 12hr at 
15 GHz, 22 GHz and 43 GHz performed on 27 February 1999, 
1 September 1999 and 20 May 2000. The 15 GHz image from 
February 1999 was used for spectral index and RM analysis by 
Gabuzda, Gomez & Agudo ( I2001t . 

The total-integration times on 0735+178 were about 
40 min at 5 GHz, 60 min at 8.4 GHz, 75 min at 15 GHz, 95 min 
at 22 GHz and 230 min at 43 GHz. 

Both left- and right-circular polarization data were 
recorded at each of the VLBA telescopes at 128 Mbps, with 
1 bit sampling and eight IFs each of 8 MHz bandwidth. The 
data were correlated at the VLBA correlator in Socorro, New 
Mexico. The calibration of the data was performed with the 
AIPS package (Fomalont 1981 1 following the standard pro- 
cedure for polarimetric observations (e.g. Leppanen, Zensus 
& Diamond 1995 1. Opacity corrections for the 22 GHz and 
43 GHz data were performed using the variation of the system 
temperature with elevation and solving for the receiver temper- 
ature and zenith opacity at each of the VLBA antennas. Opacity 
corrections for the 15 GHz data were negligible. The instru- 
mental polarization terms were determined using the parallac- 



Table 1. Observing epochs and frequencies. 



Epoch 


Date 


Freq. (GHz) 


1996.23 


26 March 1996 


8.4, 22 


1996.75 


2 October 1996 


8.4, 22 


1997.26 


6 April 1997 


5, 8.4, 15, 22 


1997.80 


18 October 1997 


8.4, 22 


1999.16 


27 February 1999 


15, 22, 43 


1999.67 


1 September 1999 


15, 22, 43 


2000.39 


20 May 2000 


15, 22, 43 



1 VLBI Space Observatory Program. 

2 The NRAO of the USA is operated by Associated Universities, 
Inc., under co-operative agreement with the National Science 
Foundation. 



tic angle dependence to separate the source polarization and the 
instrumental polarization through the feed-solution algorithm 
developed by Leppanen, Zensus & Diamond ( I1995I I. 

After the initial reduction, the data were edited, self- 
calibrated, imaged and deconvolved using CLEAN both in total 
intensity (/) and polarization (P) using a combination of AIPS 
and the Caltech DIFMAP software (Pearson et al. 19941. 

Calibration of the absolute orientation of the polariza- 
tion position angles (x) was performed by comparing the x 
values corresponding to the total polarization in the VLBA 
images for 0735+178 and various calibrators and other tar- 
get sources (PKS 0420-014, PKS 0528+134, OJ287, 3C279, 
3C 454.3 and BLLac) with those measured with the VLA 
at contemporaneous epochs. We also used the data from the 
VLA/VLBA Polarization Calibration service 3 and data from 
the 4.8 GHz. 8 GHz and 14.5 GHz UMRAO 4 monitoring pro- 
gram as a cross-check. The typical uncertainties in the x val- 
ues range between 3° and 5° at all frequencies below 43 GHz; 
at the single epoch for which we detected polarized emission in 
0735+178 at 43 GHz (May 2000), the uncertainty was larger, 
-7°. 



3. Results 

The resulting / and P images at 5 GHz, 8.4 GHz, 15 GHz, 
22 GHz and 43 GHz are presented in Figs. 1 1151 The / im- 
ages were obtained using uniformly weighted (M,v)-coverage, 
whereas the P (Stokes Q and U) images were obtained using 
natural weighting. 

The frequency coverage of our observations allowed us to 
study each region of the jet with adequate resolution within the 
inner 10 mas. The bends in the inner jet are not clearly visible 
at 5 GHz (Fig. ^ due to the relatively low angular resolution 
and the weakness of the emission near the first sharp bend at 
this frequency (due to low-frequency absorption, as reported 
by Gabuzda, Gomez & Agudo 2001). The VLBI structure at 
5 GHz is dominated by a compact core plus a broad jet ex- 
tended toward the north-east up to ~ 20 mas to ~ 30 mas from 
the core. At the other frequencies, the two previously reported 
sharp bends of ~ 90° can be seen within the inner 3 mas of the 
jet (Figs. l2l3l andl4li. Comparison of our new 8.4 GHz, 15 GHz 
and 22 GHz images and previous / images (Kellermann et al. 

3 http://www.aoc.nrao.edu/~smyers/calibration/ 

4 The University of Michigan Radio Astronomy Observatory is sup- 
ported by funds provided by the University of Michigan. 
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0735+178 at 5 GHz. April 1997. 
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Fig. 1. 5 GHz VLBA image of 0735+178 on 6 April 1997. The 
contours represent the observed total intensity, the grey scale 
the polarized intensity and the superimposed sticks show the 
orientation of the polarization electric vectors. The positions of 
the fitted Gaussian components are indicated by the crosses, 
whereas the circles (of radius equal to the FWHM of each 
Gaussian) symbolize their size. 



fT998l Gomez et al. [T9991 I200T1 Gabuzda, Gomez & Agudo 
2001} shows that the source has maintained its double bend 
structure from 1996 to 2000. However, in the last three epochs, 
the second bend seems to decrease its curvature slightly (see 
Fig.0. 

The 5 GHz and 8.4 GHz maps (Figs.[2and[2) show the po- 
larization structure extending to ~ 5 mas from the core. As was 
also reported from previous VLBI observations, the x values 
outside the first bend are aligned perpendicularly to the local 
jet axis (see e.g. Figs. |2] and [3J. The region of the first bend, 
and inwards, displays a high degree of variation in polarization 
between our observing epochs (Figs.[2]and[3j. 

3.1. Jet structure modelling 

To describe the total intensity distribution in the 0735+178 jet, 
we derived circular Gaussian model fits to all the data sets, in 
the visibility plane, using the DIFMAP software. The resulting 
parameters of the model fits, as well as the electric-vector po- 
sition angle x an d degree of polarization (m), computed as the 
X and m averages in boxes of 3 x 3 pixels around the position 
of each component, are listed in Tables[2]to[6] 

Note that a slightly different 15 GHz model fit than that 
presented by Gabuzda, Gomez & Agudo (2001 1 is presented 
here for February 1999 (Table|4j. The new model contains an 
extra component in the region between 3 mas and 4 mas from 
the core. The reason for changing this model was to make it 
more consistent with the modelling of this region at 15 GHz 
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Fig. 2. 8.4GHz VLBA images of 0735+178. Observing epochs 
are given for each image. The total intensity is plotted as con- 
tours at seven logarithmically equispaced levels, starting at the 
minimum (indicated to the right of each image) up to the 90 % 
of the corresponding / peak (also given for each map). The 
colour scale shows the distribution of the polarized intensity 
and the superimposed sticks show the orientation of the po- 
larization electric vectors. From top to bottom, the convolv- 
ing beams (shown as filled ellipses) are (1.46 x 0.90) mas 2 , 
(1.35 x 0.80)mas 2 , (1.44 x 0.79)mas 2 and (1.33 x 0.74)mas 2 
with major-axis position angles at -16.6°, -5.4°, 0.6° and 
-6.5°, respectively. The positions of the fitted Gaussian com- 
ponents are indicated by the crosses, whereas the circles (of 
radius equal to the FWHM of each Gaussian) symbolize their 



for the remaining observing epochs by considering the same 
number of components across epochs. Note, however, that this 
does not affect in any way the validity of the results presented 



by Gabuzda, Gomez & Agudo (2001 1. 

To obtain a good representation of the jet structure at 
22 GHz (epochs March 1996 and October 1996) and 43 GHz 
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Fig. 3. Same as Fig. |5] for the 15 GHz images. The location 
of the first and second bends defined in § [J are indicated 
in the image taken on April 1997. From top to bottom, the 
convolving beams are (0.85 x 0.47) mas 2 , (0.80 x 0.40) mas 2 , 
(0.73 x 0.49) mas 2 and (0.76 x 0.41) mas 2 with major axis po- 
sition angles at 0.5°, -4.7°, 1.1° and -4.8°, respectively. 
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Fig. 4. Same as Fig. |2] but for the 22 GHz images. From 
top to bottom, the convolving beams are (0.57 x 0.34) mas 2 , 
(0.57 x 0.34) mas 2 , (0.58 x 0.32) mas 2 , (0.54 x 0.30) mas 2 , 
(0.54 x 0.28) mas 2 , (0.49 x 0.34) mas 2 and (0.51 x 0.28) mas 2 
with major axis position angles at -9.5°, -9.6°, -3.0°, -6.0°, 
-4.7°, 5.4° and -2.0°, respectively. 



(epochs February 1999, September 1999 and May 2000) it was 
necessary to split the emission regions corresponding to com- 
ponents i8 and i9 (with i=C, X, U; see Figs. [2 to [5] and Tables 
|5]to|S) into several sub-components, labelled K8j, K82, Q81, 
Q82, Q83, Q9i and Q92- For a better characterization of all the 
jet components it was also necessary to fit the extended emis- 
sion observed in the outer regions of the jet (modelled as com- 
ponents iE with i=C, X, U, K; see Figs. [2 to 0] and Tables |2] 
to |3- Note that these components should not be taken as rep- 
resenting real distinct features in the jet, but instead extended 
regions with weak flux densities. 



The errors in the DIFMAP model-fit parameters (total flux 
density, position and size of the components) were obtained 
using the DIFWRAP software, developed by Lovell (2000 1. 
The errors were determined by perturbing the fit parameters 
until the resulting residual maps were unacceptable. These per- 
turbed quantities defined the upper and lower limits for the er- 
rors of the fit parameters. To take into account the interrela- 
tion between the fit parameters, we also allowed for simultane- 
ous variations of all four parameters for each component when 
determining the final uncertainties. This resulted in somewhat 
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0735+1 78 at 43GHz 
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Fig. 5. Same as Fig. |2]but for the 43 GHz images. From top to 
bottom, the convolving beams are (0.31 x 0.16) mas 2 , (0.29 x 
0.20) mas 2 and (0.28 x 0.15) mas 2 with major axis position an- 
gles at -12.8°, 7.8° and -3.4°, respectively. 

larger but more conservative error estimates (see Lovell |2000) 
and http://www.vsop.isas.ac.jp/survey/difwrapi. 

The;f errors were computed from the uncertainties in the^ 
calibration for each map added in quadrature to the rms noise 
within the 3x3 pixel box in which x was averaged. The m 
uncertainties were obtained as the propagated errors from the 
/ and P measurements within the averaging box, for which the 
rms noise of / and P was taken as uncertainties of their corre- 
sponding measurements. 

4. Discussion 

4.1. Kinematics of model components 

Figure ^represents the positions of all the model components 
in Tables |2| to |6| within the inner 6 mas of the jet for which 
reliable position estimates were obtained. This figure shows a 
clustering of components around some jet locations, marked by 
dashed circles. This clustering was rather stable during a period 
of four years and is very unlikely to be produced by chance. 
This kind of component distribution was also shown by Gomez 
et al. (2001 1 for a different series of data spanning from 1995 to 



1998. The clustering implies a much slower scenario than re- 
ported for epochs before 1993 (Gabuzda et al. 1994 and refer- 
ences therein), with proper motions in the range 0.28 mas yr~' 
(7 c) to 0.48 mas yr~' (12 c). Such proper motions would be ob- 
served as systematic shifts of the positions of the model com- 
ponents by 1.1 2 mas to 1.92 mas during the time range of our 
monitoring (~ 4yr), which are not observed in any of the jet 
regions. 

Figure presents the distance from the core as a function 
of time for the innermost (and most reliable) 8.4 GHz, 15 GHz, 
22 GHz and 43 GHz model components in Tables [5] to |6] and 
those published by Gomez et al. (2001 1. The 5 GHz data in 
Table|2]are not plotted due to the much lower resolution at this 
frequency, which can produce confusion in the identification of 
components. The evidence of slow proper motions along the jet 
allowed us to ensure the most reliable identification of compo- 
nents as the one which relates a jet feature with those nearest 
in position at different observing frequencies and at adjacent 
observing epochs. 

It is worth to note the significant discrepancy between the 
identification of components presented here and the one per- 
formed by Kellermann et al. (2004) through their quasi-annual 
15 GHz observations of 0735+178 from 1995 to 2000. This 
discrepancy is most probably caused by the poorer time and 
frequency sampling of Kellermann et al. and also by the differ- 
ent model-fitting method used by them, which only considered 
a maximum of three jet components per epoch and was per- 
formed in the image plane. Figure shows that a set of only 
three jet components does not suffice to optimally fit the whole 
VLBI structure of 0735+178 at 15 GHz. 

Our fits of the radial positions of components i8, i9 and 
ilO show slow radial motions outward from the core with 
speeds of (0.023 + 0.004) mas yr 1 , (0.048 + 0.004)masyr~ 1 
and (0.026 + O^^masyr- 1 ([0.58 + 0.10] c, [1.21 ± 0.10] c 
and [0.65 + 0.53] c), respectively. The proper motions of i8 
and i9 were interpreted as consistent with a quasi-stationary 
behaviour by Gomez et al. (2001 1 for a shorter time range 
(from 1995 to 1998) than the one covered by our observa- 
tions. Figure [6] shows that these two model components are 
more consistently fitted as moving features if our last three 
observing epochs (from 1999 to 2000) are taken into account. 
This implies that both i8 and i9 slightly speeded up after 1998, 
which is more evident for i9 (Fig. |6j. It should also be noted 
that our estimates provide only the speeds of the average po- 
sition of the components during our observations. The unac- 
celerated proper motion of 0.14 + 0.03 mas yr -1 estimated by 
Homan et al. i 



(2001 1 for i8, through their bimonthly monitoring 
program during 1996, accounts for the irregular motion of this 
component around its average position. This rapid, but short 
distance, motion explains the clumping of components around 
their mean positions, which are most likely produced by the 
changes of their brightness distributions. 

The fits for components il, i3 and i4 show slight upstream 
radial motions. Although these motions could, in principle, 
happen in actual relativistic jets under certain conditions (see 
e.g. Agudo et al. 2001 1, we prefer to consider these three 
components to be quasi-stationary in view of the large un- 
certainties in their positions and apparent motions (with val- 
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Table 2. 5 GHz model fits for 6 April 1997. Table 3. 8.4 GHz model fits. 



Comp. 



/ 

(mJy) 



(mas) 



e 
O 



FWHM 
(mas) 



X 
(') 



C8 

C5+6 

C3+4 

CI 

CE1 

CE2 



249±12 
455±8 
100±13 
47±11 

19±9 
70±40 

-81 



0.00 
0.78±0.09 
1.90±0.15 
3.88±0.15 
5.38±0.60 
10.33±1.19 
-21.79 




67±30 
41±5 
65±6 
66±10 
73±6 
-76 



<0.43 
0.82±0.08 
1.55±0.23 
1.18±0.51 
1.59±0.93 
6.06±1.67 
-16.40 



-23±3 2.2±0.3 

-14±3 4.1±0.4 

-1±3 10±2 

-11±4 8±5 



ues of [-0.11 + 0.19] mas yr 1 , [-0.05 + 0.09] mas yr" 1 and 
[-0.07 + 0.07] mas yr~', respectively), which are in fact con- 
sistent with no motion. Further evidence for considering these 
components as quasi-stationary is based on their stable struc- 
tural position angles, with average values of 66° + 7°, 65° + 5° 
and 64° + 5°, respectively, and their stable total flux densities 
(see §0. 

The behaviour of the components related to the second 
bend, i5 and i6, is very different. Inspection of Fig.[3]reveals 
a smooth change in the region of the second bend, which pro- 
gressively softens from 1997 to 2000. Figure [6] also shows a 
significant dispersion for the position of component i6 (from 
March 1996 to October 1997), which is larger than its typical 
position uncertainty (especially at 8.4 GHz). Fig. [8]represents 
the structural position angles of components i5, X5 1 and i6 dur- 
ing the time covered by our observations. We can see that i6 
moved from a structural position angle of ~ 38° at the begin- 
ning of 1996 to ~ 20° from the end of 1996 to the beginning 
of 1997. During 1999, it returned to ~ 30°, then finally became 
undetectable in our last observing epoch (May 2000). Although 
the variations are less significant, the structural position angle 
of i5 also changed from ~ 47° at the beginning of 1996, to 
~ 52°, during 1997, to ~ 37° at the beginning of 1999, and 
then back to almost the initial position angle at the beginning 
of 2000. During the time period of maximum angular separa- 
tion between i5 and i6 (during 1997) a new component (X5i) 
was fitted at intermediate angular positions of 35° + 2° (in April 
1997) and 32° + 2° (in October 1997). Although it was neces- 
sary to introduce component X5j for these two epochs in order 
to obtain a good model fit, it seems very unlikely that this is 
related to the appearance of a new real feature in the jet. We 
consider X5i to represent emission from the underlying jet in 
the gap region produced as a consequence of the angular sepa- 
ration of i5 andi6. 

The behaviour of i5 and i6 between 1996 and 2000 is re- 
lated to a "softening" of the second bend, which is consis- 
tent with the fairly rectilinear east-north-east jet structure re- 



ported by the 15 GHz image of Lister & Homan (2005) taken 
on November 2002. This softening is produced by an apparent 
backward motion of component i6. Again, this should not be 
understood as a physical motion of jet plasma back to the core, 
but rather as induced by a decrease in the total flux density of 
the bent area, shifting the peak of the emission (associated with 
the component) to positions closer to the core until component 
i6 became undetectable by May 2000 (see Tables |2] to [6] and 

§E3. 



Comp. 



/ 

(mJy) 



(mas) 



9 
C) 



FWHM 
(mas) 



X 

C) 



26 March 1996 



Core 
X8 

X6 
X5 
X4 
X3 
XI 
XE 



590±39 
629±24 
161±24 

16±9 
23±12 
21±15 
20±15 

-77 



0.00 
0.76±0.01 
1.45±0.02 
2.68±0.27 
3.51±0.24 
4.41±0.47 
5.43±0.66 
-11.08 




77±2 
38±2 
47±5 
65±9 
65±8 
66±13 
-72 



0.20±0. 1 1 
0.42±0.06 
0.91 ±0.22 
<0.54 
<0.61 
<1.21 
<2.10 
-8.06 



-0.8±3 7±1 

5±3 23±11 

-16±3 14±8 



2 October 1996 



Core 

X8 
X6 
X5 
X4 
X3 
XI 
XE 



396±4 
434±6 
110±11 
36±8 
15±8 
18±8 
19±12 
-62 



0.00 
0.80*0.01 

1.39±0.02 
2.24±0.10 
3.39±0.24 
4.11±0.16 
5.32±0.35 
-10.8 




77±1 
27±1 
46±2 
64±7 
64±4 
66±6 
-73 



0.25±0.04 
0.53±0.02 
0.86±0.03 
0.99±0.15 
1.00±0.43 
0.72±0.43 
1.26±0.53 
-7.09 



-47±3 
-49±3 
-15*3 
13±3 
-21±3 
-36±3 



2.4±0.2 

2.7±0.3 

8.6±0.8 

25±5 

19±6 

17±6 



6 April 1997 



Core 

X8 
X6 
X5, 
X5, 
X4 
X3 
XI 
XE 



350±4 
310±6 
52±8 
31±7 
19±5 
23±7 
12±4 
17±8 
-56 



0.00 
0.85±0.01 
1.46*0.05 

2.12±0.08 
2.52±0.09 
3.48±0.12 
4.17±0.12 
5.32*0.25 
-11.04 




75*1 
17*2 
35*2 
52*3 
64*3 
65*2 
67*5 
-73 



0.22*0.03 
0.67*0.02 
0.60*0.10 
0.67*0.15 
0.53*0.34 
1.08*0.26 

<0.56 
1.39*0.73 

-6.41 



-3*4 
-25*3 
-11*3 
1*3 
2*3 
-11*3 
-23*3 



1.1*0.3 
3.9*0.6 

15.8*1.7 
21*4 
19*6 
17*7 
14*8 



18 October 1997 



Core 
X8 
X6 
X5i 

X5> 
X4 
X3 
XI 
XE 



391*2 
222*2 
40*4 
28*4 
37*4 
18*7 
18*6 
10*9 
-67 



0.00 
0.83*0.01 
1.61*0.05 
2.11*0.10 
2.57*0.09 
3.53*0.17 
4.33*0.17 
5.72*0.67 
-11.28 




73*1 
18*1 
32*2 
52*2 
67*4 
65*3 
66*9 
-73 



0.18*0.01 
0.71*0.03 
0.51*0.18 
0.48*0.23 
0.78*0.24 
0.76*0.44 
0.75*0.44 
1.37*0.64 
-6.98 



80*3 
-84*3 
-52*3 
-34*3 
-20*3 
-44*3 
-47*3 



3.2*0.2 

4.9*0.5 

14.5*1.4 

24.3*2.2 

23*3 

14*5 

16*5 



4.2. Flux density evolution of model components 

Fig. 13 presents the 5 GHz, 8.4GHz, 15 GHz and 22 GHz to- 
tal flux density evolution for the core and components ilO, 
i9, i8, i6 and i5 during the time span covered by our obser- 
vations. Components not included in this figure, X5i; i4; i3 
and il, showed little 8.4 GHz total flux density evolution with 
weighted means of (29 ± 4)mJy, (17 + 4) mJy ([17 + 8] mJy at 
15GHz),(16 + 2)mJy([9 + 5]mJyatl5GHz)and(ll+3)mJy, 
respectively. 

This figure shows that the total flux density of i6 at 8.4 GHz 
decreased from ~ 160mJy at the beginning of 1996 to ~ 
10 mJy at the end of 1999, while the nearby feature i5 displayed 
a quite stable total flux density of (35 + 6) mJy at 8.4 GHz and 
(20 + 5)mJy at 15 GHz; i5 was still detected in May 2000, 
whereas i6 had become undetectable by that time. 

A similar trend as for i6 was observed for component i8 
at the first bend, which quasi-monotonically decreased its to- 
tal flux density by ~ 500 mJy at 8.4 GHz during the same pe- 
riod. Simultaneously, the total flux density of the nearby com- 
ponent i9 remained almost constant at ~ 100 mJy from 1996 
to 1999. Component ilO likewise maintained a very stable total 
flux density of (74 + 7) mJy at 43 GHz from February 1999 to 
May 2000. 
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Fig. 6. Component positions and corresponding uncertainties (symbolized by the size of the crosses) for the model fits presented 
in Tables |2 to |6] The dashed circles mark the various component clumps observed at each frequency. Components with large 
uncertainties in the fit parameters (for which no errors are given in Tables[2]to[6) and those with an uncertain identification (see 
Tables0]and[5} are not plotted. 
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Fig. 7. Distance to the core as a function of time for the inner- 
most (within < 7 mas) model components at 8.4 GHz, 15 GHz, 
22 GHz and 43 GHz. Note that, for the sake of clarity, only 
components observed during more than two epochs have been 
plotted. The epochs used are those in Tables 0to [6] as well as 
August 1996, January 1997 and May 1998 from Gomez et al. 
(2001). Linear fits of the radial positions for some model com- 
ponents are also shown (solid lines). The dashed lines connect 
the radial positions of components i5, X5j and i6. Components 
with large uncertainties in the fit parameters (for which no er- 
rors are given in Tables |3 to [6} and those with an uncertain 
identification (see Tables |4]and[5j are not plotted. 
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Fig. 8. Structural position angles as a function of time for com- 
ponents i5, X5i and i6. The points shown are the same as those 
in Fig. 



While the total flux density evolution of i8, i9 and ilO was 
monotonic for most of the time between 1996 and 2000, the 
core displayed significant rapid fluctuations on typical time- 
scales of several weeks and with amplitudes of ~ 20 % around 
400 mJy. 

Note that, whereas the total flux density of the core fluctu- 
ated and the flux densities of i9 and ilO remained almost con- 
stant, the behaviour of i8 coincided with the single-dish inte- 
grated total flux density evolution of 0735+178 between early 
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Table 4. 15 GHz model fits. 



Table 5. 22 GHz model fits. 
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5 April 1997 












26 March 1996 






Core 


292±7 


0.00 





0.15+0.02 


-9±4 


0.6±0.5 


Core 


353±9 


0.00 





<0.10 






U9 1 


15±8 


0.40±0.23 


53±23 


<0.38 






K9 


73±9 


0.20±0.03 


76±18 


<0.32 






US 


236±18 


0.90±0.06 


75±2 


0.67±0.03 


-39±3 


2.6±1.3 


K8i 


232±13 


0.75±0.01 


88±1 


0.34±0.04 


60±4 


4+3 


U6 


47 ±9 


1.50±0.10 


19±2 


0.69±0.18 


-25±3 


18±7 


K7 


122±16 


0.95±0.05 


53±4 


0.68±0.10 






U5 


24±12 


2.25±0.14 


43±4 


0.76±0.34 


-1±3 


22±15 


K8, 


93±13 


1.08*0.04 


70±3 


0.34±0.07 






U4 1 


10±10 


2.92±0.33 


61±8 


0.53±0.39 






K6 


-5 


-1.57 


-13 


-0.01 






U3 


11±9 


3.96±0.31 


66±8 


±1.15 






K5 1 


45±34 


1.90±0.20 


35±8 


1.08±0.34 






UE 


-14 


-4.83 


-64 


-1.45 






K4 


-8 


-3.12 


-61 


-0.68 










27 


February 1999 






K3 


-10 


-4.07 


-65 


-0.67 
-0.90 






Core 


417±2 


0.00 





0.09±0.02 


-69±3 


1.3±0.2 














U9 


83±4 


26±0 01 


94±5 


22±0 13 


-68±3 


2 0±0 4 












U8 


47±9 


0.96±0.05 


80±5 


0.55±0.18 


-46±3 


11±4 


Core 


217±7 


0.00 





<0.10 


14±3 


2.2±1.2 


U6 


25±15 


1.47±0.25 


29±9 


0.91±0.32 






K9 


33±8 


0.32±0.05 


88±12 


0.26±0.05 






U5 


13±9 


2.32±0.25 


37±6 


0.60±0.44 


-11±3 


36±17 


K8, 


70±10 


0.86*0.02 


91±2 


0.30±0.09 


-72±3 


6+4 


U4 


18±13 


3.02±0.35 


58±9 


1.08±0.66 






K8 2 


153±10 


0.91 ±0.02 


67±3 


0.59±0.08 






U3 


13±13 


3.83±0.44 


66±8 


1.03±1.01 






K6 


26±14 


1.51±0.16 


22±6 


0.68±0.31 


-9*3 


34±30 


UE 


~6 


-5.31 


-65 


-1.26 






K5 


-22 


-2.31 


-47 


-1.28 










1 September 1999 












5 April 1997 








Core 


365±4 


0.00 





0.12±0.03 


-38±5 


3.6±0.3 


Core 


184±11 


0.00 





<0.18 


-6±4 


2.7±1.4 


U9 


176±5 


0.32±O.Ol 


76±3 


0.24±0.07 


-39±5 


7.0±0.7 


K9 


47±11 


0.20±0.05 


67±45 


<0.35 


5*3 


4.9*2.4 


US 


38±8 


0.88±0.06 


79±6 


0.57±0.16 


-44±5 


13±6 


K8 


171±28 


0.92±0.03 


77±2 


0.69±0.06 






U6 


9±8 


1.43±0.39 


28±9 


0.48±0.44 






K6 


43±24 


1.50±0.17 


23±8 


0.79±0.20 


4*3 


36±29 


U5 


24±9 


2.18±0.23 


41±6 


0.89±0.28 


-4±5 


21±18 


K5 


-19 


-2.41 


-46 


-0.93 


—2 


-86 


U4 


17±10 


3.33±0.39 


61±6 


0.92±0.44 






KE[ 


-18 


-4.08 


-64 


-1.09 






U3 

UE 


8 ±6 
-6 


3.95±0.39 
-5.04 


66±6 
-67 


<1.31 
-1.20 










18 October 1997 






Core 


244±6 


0.00 





<0.13 


18*3 


5.2±0.5 






2u May ZUUl 








K9 


43±8 


0.19±0.03 


85*17 


<0.21 


38±4 


6.0*1.1 


Core 


523±6 


0.00 





0.12*0.02 


-60±5 


2.0±0.2 


K8 


107±16 


0.89*0.06 


76*5 


0.69±0.05 






U9 


144±9 


0.33±O.Ol 


79±7 


0.30±0.05 


-63±5 


2.9±0.6 


K6 


34±20 


1.52±0.23 


23±11 


0.79*0.35 






US 


86±8 


0.84±0.02 


76±3 


0.54±0.05 


-70±5 


10±3 


K5 


-25 


-2.36 


-42 


-1.08 






U5 


30±17 


2.13±0.34 


43±15 


1.31±0.91 






KE[ 


-20 


-3.87 


-65 


-1.37 






U4 
U3 

UE 


~7 
-15 
-3 


-3.16 
-3.97 

-5.72 


-62 
-66 
-69 


-0.57 
-0.88 
-1.03 










27 


February 1999 






Core 
K9 


414±9 
76±18 


0.00 
0.24±0.03 



85±23 


<0.19 
0.25*0.11 


-42±3 


0.8±0.4 
















K8 


33±28 


0.95±0.14 


85±10 


0.49±0.24 






'Uncertain identification of this model 


component 






KE 2 


-60 


-1.80 


-42 


-1.73 
























1 September 1999 








Core 


345±11 


0.00 





0.09±0.02 


-32±3 


2.1*0.6 
















K9 


145±14 


0.32±0.02 


75*11 


0.26±0.10 


-29±3 


7.3*1.7 
















K8 


24±18 


0.87±0.15 


83±14 


0.51*0.26 




















KE 2 


-47 


-2.31 


-49 


-2.06 






1996 anc 


1 late 1997 at 4.8 GHz, 8 GHz and 14.5 GHz (Aller et 
and at 22 GHz and 37 GHz (Fig.^|and Terasranta et 
. The data from both single dish monitoring programs 






20 May 2000 








al. [19991 
al.l2004l 


Core 
K9 

K8 


520±7 
134±19 
92±16 


0.00 
0.29±0.02 
0.76±0.07 



81±5 
75*6 


0.12±0.01 
0.31*0.07 
0.61*0.09 


-52±3 
-78±3 


1.7*0.3 

3.3±1.6 


showed c 


monotonic decrease in total flux density of 


~500mJy 


KE 2 


-42 


-2.91 


-57 


-2.27 







at all frequencies during this time range, which is precisely the 
amount of flux density decrease observed for i8. 

This, together with the fact that i8, i9, ilO and the core are 
typically responsible for more than the 90 % of the total emis- 
sion of the source at high frequencies (22 GHz and 43 GHz) and 
more than 80 % at lower frequencies, shows that the decrease in 
total flux density measured by the single-dish monitoring was 
determined primarily by the flux density decrease of i8. 

4.3. Doppler factor and geometry of the inner jet 

The highly bent VLBI structure of 0735+178 suggests that the 
source is oriented at a small angle to our line of sight. This is 
also supported by the large superluminal motions measured in 
the jet before 1993, and by the observed high degree of vari- 
ability for the core (see 34.21 

It is possible to constrain this angle using our proper motion 
measurements and previous estimates of the jet's Doppler fac- 



1 Uncertain identification of this model component. 



tor 5 (6). Giiijosa & Daly ( 1996 1 computed a 5 eq =» 7 (estimation 
of the Doppler factor under the assumption that the magnetic 
field and the relativistic particles of the plasma are in equipar- 
tition). Using different assumptions, Lahteenmaki & Valtaoja 
(1999) reported 6 mr ~ 4 (Doppler factor estimated from the 
variability of the long term radio light curves). Other estimates 
of 6 for 0735+178 in the literature lie between the previous val- 
ues of 5 eq and 5 var . To test these values, we have computed the 
variability Doppler factor 8 mr of the core in 0735+178 through 
the expression defined by Jorstad et al. J2005l > and using our 
own data. The result gives 6 var ~ 5, which agrees well with the 
previous estimates. 



5 S = [T(l -pcos6)Y l where T = (1 - /J 2 )"" 2 is the Lorentz factor, 
P is the speed in units of the speed of light and 6 is the angle between 
the direction of the flow and the line of sight. 
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Table 6. 43 GHz model fits. 
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27 February 1999 



Core 


315*12 


0.00 





0.07±0.01 


Q10 


68±9 


0.08±0.03 


123±13 


<0.22 


Q9, 


62±27 


0.21±0.02 


64±14 


0.21±0.04 


09- 


~7 


-0.40 


-98 


-0.02 


08, 


-5 


-0.62 


-92 


-0.10 


08, 


-11 


-0.79 


-75 


-0.28 


Q8 3 


~8 


-1.08 


-86 


-0.22 



X 

O 



1 September 1999 


Core 

Q10 
Q9, 
Q9 2 


257±14 
79±18 
43 ±29 
59±43 


0.00 
0.12±0.02 
0.26±0.06 
0.42±0.10 


<0.08 
70±21 <0.20 
58±14 <0.34 
83±19 0.29±0.22 


20 May 2000 



252±12 
163 ±9 
82±12 

162±34 

-14 
-18 



0.00 

0.06±0.01 

0.13±0.01 

0.33±O.O8 

-0.76 

-1.10 




103±25 
84±11 

77±16 
-71 
-81 



0.07±0.04 
<0.08 
<0.21 

0.47±0.07 
-0.20 
-0.32 



-51 ±7 
-59±8 

-74±8 



1.4*0.5 
1.4*0.6 
3.1*1.1 
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Fig. 10. Integrated 22 GHz and 37 GHz total flux density light 
curves of 0735+178 from 1996 to 2000. Data obtained by 
the monitoring program of extragalactic radio sources of the 
Metsahovi radio observatory (Terasranta et al. 2004 1. 
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Fig. 9. 5 GHz, 8.4 GHz, 15 GHz and 22 GHz and 43 GHz total 
flux density evolution for the Core, ilO, i9, i8, i6 and i5 compo- 
nents. The epochs shown are those in Tabled as well as August 
1996, January 1997 and May 1998 from Gomez et al. (2001). 
Each frequency is represented by a different symbol, as indi- 
cated in the right top inset. Note the different total flux density 
scale used for the plots of i6 and i5. 



We then assume that the true S for the innermost jet in 
0735+178 ranges between 4 and 7, and use our observed 
proper motions of components i8, i9 and ilO (which essen- 
tially range between 0.1 c and 1.3 c) to estimate the possible 
ranges of viewing angles (6) and plasma bulk Lorentz fac- 



tors (T) for the first mas of the jet, using the relation /3 ap = 
(J3sinff)/{1 - ficosff) where /3 ap is the measured apparent proper 
motion (e.g., Pearson & Zensus 1987 1 together with the defini- 
tion of the Doppler factor (see above). This constrains 8 to val- 
ues between 0° and 9°, as small as expected, and constrains F to 
values between 2 and 4. Note that the speeds used for these es- 
timates are associated with jet features that are related to bends 
(at least for i8; see Figs. 0] and and we have not observed 
fast superluminal speeds such as those reported before 1993. 
Therefore, our estimates of 6 and T should be related to the 
pattern speeds of the innermost bent jet structure. 

4.4. Polarization evolution of model components 

Fig.llllrepresents the^ values at the locations of the model fit 
components in Tables [2] to [6] In agreement with Gomez et al. 
(1999), this figure shows that the polarization electric vectors 
were predominantly perpendicular to the local direction of the 
northern jet (outwards from the location of i6) at all the observ- 
ing epochs. 

For the region between the core and the first strong bend, 
X shows a clear time dependence, which is more clearly shown 
by the plots of^ as a function of time for i8 and the core (Fig. 
I12> . Note that the x behaviour of i9 is very similar to that of 
the core. This is because these two components belong to the 
same polarization region, although the model fits for the total 
intensity jet are able to separate them. While the core shows 
irregular polarization variability from 1996 to 2000 (as is ex- 
pected from its high level of total flux density variability), the 

X for i8 linearly increases at all frequencies between 120° 

in March 1996 to ~ -30° in April 1997. Furthermore, it expe- 
rienced an abrupt rotation by ~ -60° in October 1997 and then 

maintained a much more stable value of 45° during the last 

three observing epochs, with a small decrease to ~ -25° ob- 
served in May 2000. 

The ~ 90° monotonic rotation of i8, from 1996 to 1998, 
shows evidence of correlation with its total flux density evo- 
lution and also with the integrated total flux density evolution 
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Fig. 11. Time-dependent^ distribution at the location of the 
model fit components from Tables [2] to [6] The lines represent 
the orientation of x, whereas the observing frequency for each 
component is shown by a different symbol (as indicated in the 
left bottom panel). 



• 5 GHz 
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Fig. 12. x as a function of time for i8 (top) and the core (bot- 
tom). The observing frequency for each component is repre- 
sented by a different symbol, as indicated in the top plot. 



of 0735+178 (see 34. 2> during the same time range. This rota- 
tion can not be produced by time dependent Faraday rotation, 
which would progressively change the relative ;^ orientation of 
i8 for the different frequencies. A possible explanation for this 
behaviour could be either the passage of a shock through the 
first bend or the interaction of this bend with a clumpy external 
medium, as proposed by Gabuzda, Gomez & Agudo (2001 1, 
see also § [l] We favour the jet-external medium interaction 
since, although i8 moved outward from the core, it did so with 
a slow subluminal speed not typical of moving shocks in rela- 
tivistic jets and it displayed a stable position and seemed to be 
well "attached" to the first bend of the underlying jet. In addi- 
tion, a shock passing through the first bend would be expected 
to produce a clockwise rotation of x (f° r a plane shock per- 
pendicular to the jet axis and moving following the jet bending 
towards the north), contrary to what was shown by our obser- 
vations. 

Note also that the behaviour of i8 is quite different from 
that of i6, which maintained a stable x orientation during the 
time covered by our observations, although its total flux density 
monotonically decreased (as for i8). 



5. Summary and conclusions 

We have investigated the kinematic, total and polarized flux 
density evolution of roughly the inner 6 mas of the VLBI jet in 
0735+178 from 1996 to 2000. 

Components i 1 , 13 and i4 displayed essentially stationary 
positions during the time spanned by our observations. They 
also have stable / and P flux densities and^ values. This sug- 
gests that the region in which these components are located is 
mainly composed of the quiescent underlying jet, whose dy- 
namical time-scales seem to be longer than those covered by 
our observations. 

For the innermost jet, we have found subluminal or slow 
superluminal proper motions with apparent speeds of (0.58 + 
0.10) c, (1.21 +0.10)c and (0.65 +0.53)c for components i8, i9 
and ilO, respectively. The variability of the total flux density of 
the core allowed us to estimate its Doppler factor 5 ~ 5. Using 
these measurements and previous estimates of the jet Doppler 
factor, we have computed the possible values for the Lorentz 
factor T, which ranged from 2 to 4, and the angle to the line 
of sight 9 to be < 9°. As we have not detected fast superlumi- 
nal motions, such as those reported for this source before 1993, 
and the apparent speeds used to estimate 9 and F were based on 
the slow proper motions of jet features near the jet bends, our 
estimates of F should be related to the typical pattern speeds 
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of the innermost bent jet structure. Note that Lorentz factors in 
bent regions similar to those estimated for 0735+178 have been 
shown by three-dimensional numerical simulations of relativis- 
tic helical jets (Aloy et al. 2003 1. 

Therefore, with the exception of i6 (whose position 
changed significantly, moving first north and then back south, 
leading to a decrease in the curvature of the second bend of the 
jet), all the jet components traced during the time range cov- 
ered by our monitoring were either quasi-stationary or showed 
slow motions, contrary to the situation observed before 1993 
when 1 c < /3 ap < 12 c (Gabuzda et al. ( 1994i and references 
therein). 

A possible explanation for this discrepancy was suggested 
by Gomez et al. (2001 1, who proposed a scenario in which the 
jet experienced different regimes before 1993 and after 1995. 
Before 1993, rapid superluminal motions (7 c < f3 ap < 12 c) 
were reported, the structure of the jet remained mainly recti- 
linear (Gabuzda et al. 119941 and Gomez et al. 2001 1 and the 
radio flux densities from the UMRAO and Metsahovi monitor- 
ing programs were the highest ever detected in 0735+178 (up 
to ~ 5 Jy, see Aller et al. [T9^9l and Terasranta et al. 1277041 . This 
could be called an "active" regime of the jet in 0735+178. From 
1995 to 2000, in the "quiescent" regime, only slow motions 
were detected (see 34. 1> . the source showed a highly curved 
structure with two sharp bends and its integrated flux densities 
between 5 GHz and 37 GHz decreased roughly linearly up to 
~ 1997 to ~ 1998, reaching its lowest known values ( < 1 Jy, 
Aller et al. lT999l and Terasranta et al. l2004l 

This scenario is supported by the 15 GHz VLBA images 
of 0735+178 obtained by the MOJAVE program 6 for epochs 
after November 2002 (Lister & Homan 2005 1. In agreement 
with the softening of the second bend during the time covered 
by our observations, these new images show a fairly rectilin- 
ear jet structure oriented to the north-east (as reported for early 
VLBI images before 1993). The typical double-bend structure 
displayed by 0735+178 between 1995 and 2000 is no longer 
evident. In addition, the 15 GHz integrated total flux density of 
the source in November 2002 increased by 100 % with respect 
to our previous observation at 15 GHz (in May 2000), which 
seems to be accompanied by the appearance of rapid superlu- 
minal motions within the innermost four mas of the jet. 

In addition, the optical (V, R and / band) light curves of 
0735+178 from January 1995 to March 2001 reported by Qian 
& Tao (2004) also agree with the proposed scenario. The source 
was observed to decrease its optical flux density from 1995 to 
January 1998, when it reached one of its faintest levels ever 
observed (with a magnitude in V, my - 16.68). After that, 
0735+178 became brighter again to reach a new maximum in 
February 2001 (my - 14.54). This behaviour is consistent with 
the one showed by the radio light curves in Fig. [K)] Hanski, 
Takalo & Valtaoja (2002i, through the analysis of a different 
optical data set from 1988 to 1998 together with the radio light 
curves from the Metsahovi monitoring program, showed that 
the optical flux density level was always high for 0735+178 
when the modelled radio flux density was at the peak of a flare. 



http://www.physics.purdue.edu/astro/MOJAVE/ 



Confirmation of superluminal motions from the data of the 
MOJAVE program will provide added support for the proposed 
bimodal scenario. We could therefore consider the possibility 
of a cyclic behaviour for 0735+178 in which periods of en- 
hanced activity, with ejection of superluminal components, are 
followed by epochs of low activity, with a highly twisted jet 
geometry, and vice-versa. During the high activity period the 
ejection of superluminal components, associated with strong 
shocks, would lead to an increase of the broad band emission 
of the source, enhanced by the large Doppler boosting factors 
expected for this source given its small viewing angle. An in- 
crease in the momentum of the component may lead to a more 
rectilinear, quasi-ballistic, motion of the component, and there- 
fore for the apparent jet geometry. Subsequent low activity pe- 
riods could reveal the quiescent jet, characterized by a highly 
twisted geometry. This would most probably require a jet with 
a changing direction of ejection, since otherwise we would ex- 
pect the quiescent jet to follow the straight funnel left by the 
previous period of high activity. 
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